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ABSTRACT 

We present a light curve model for the super-Eddington luminosity phase of five classical novae observed 
with IUE. Optical and UV light curves are calculated based on the optically thick wind theory with a reduced 
effective opacity for a porous atmosphere. Fitting a model light curve with the UV 1455 A light curve, we 
determine the white dwarf mass and distance to be (1.3 M©, 4.4 kpc) for V693 CrA, (1.05 M , 1.8 kpc) for 
V1974 Cyg, (0.95 M , 4.1 kpc) for V1668 Cyg, (1.0 M , 2.1 kpc) for V351 Pup, and (1.0 M Q , 4.3 kpc) for 
OS And. 

Subject headings: novae, cataclysmic variables — stars: individual (V693 Coronae Australis, VI 974 Cygni, 
V1668 Cygni, V351 Puppis, OS Andromedae) 
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1. INTRODUCTION 

The super-Eddington luminosity is one of the long stand- 
ing problems in the theoretical study of classical novae (e.g., 
Friediung 2004} for recent summary). Super-Eddington 
phases last more than a few to several days, and their peak 
luminosities often exceed the Eddington limit by a factor of a 
few to several (e.g. jDella Valle & Livioll995l and references 
therein). 

It is difficult to reproduce such a super-Eddington luminos- 
ity in evolutional calculations of nova outbursts. Dynami- 
cal calculations show that the super-Eddington phase appears 
only in a very short time, or does not appear at all, apart 
from numerical difficulties that often prevent accurate cal- 
culation of the photospheric luminosity and v i sual magni- 
tude ( e.g.JSparks et al 11 1 9781; iPrialnik et all 1978b iNariai et all 
ir^lSterrfield et al. II1985T 1986[lPolitano et allll995l) . 

Recently, IShaviv I d2001l 120021) presented an idea on the 
mechanism of the super-Eddington luminosity. Shortly af- 
ter hydrogen ignites on a white dwarf, the envelope becomes 
unstable to develop a porous structure in which the effec- 
tive opacity becomes much smaller than the normal opacity 
for uniform medium. Corresponding to the reduced effective 
opacity, the effective Eddington luminosity becomes larger. 
Therefore, the diffusive energy flux can exceed the Eddington 
value for uniform medium, even though it does not exceed the 
effective Eddington lu minosity. 

Based on this idea, iKato & Hachisu I ( 120051) presented a 
light curve model for the super-Eddington phase of classi- 
cal novae. They assumed a reduced opacity in early phases 
of nova outbursts and reproduced the optical light curve of 
V1974 Cyg. This is the first theoretical model for the super- 
Eddington light curves. In the present paper, we apply the 
same method to other classical novae to examine whether or 
not this idea is applicable to the different speed class of novae. 

The International Ultraviolet Explorer (IUE) satellite 
obser ve d a nu m ber of n o va outb ur sts (e.g., [ Cassatell a et al. I 
[19791 l200l l2004al 12001 IStickland et al. I 119811) . 
Cassa tella et al. I ((2002) presented 1455 A continuum 
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light curves for twelve novae and showed that the duration 
of the UV outburst is a good indicator of the speed class of 
novae, i.e., a faster nova shows a shorter duration of the UV 
outburst. This 1455 A light curve is an important clue in 
modeling the super-Eddington phase dKato & Hachisu |[2005l) 
and also a useful too l in estimating the white dwar f mass and 
distance to the star (lHachisu & Kato Il200l |2006|) . Here, we 
make a model of the super-Eddington phase for five Galactic 
novae, V693 CrA, VI 974 Cyg, VI 668 Cyg, V351 Pup, and 
OS And, because their 1455 A light curves are available 
from the beginning of the super-Eddington phase to the UV 
decay. 

Section [2] gives a brief description of our method for UV 
light-curve fittings. The light curve analyzes for individual 
objects are shown in $3|-[7] In $8]we summarize our results. 

2. LIGHT CURVE MODEL 
2.1. Optically Thick Wind Model 

After a thermonuclear runaway sets in on an accreting 
white dwarf (WD), the photosphere greatly expands to R p h > 
100 Rq. The optical luminosity reaches a maximum value, 
which often exceeds the Eddington limit. After that, the pho- 
tosphere moves inward whereas the envelope matter goes out- 
ward. The wind mass-loss begins in the very early phase of 
the outburst and continues until the photospheric temperature 
rises to log7J,h ^5.2-5.6. The enve lope mass decreases ow - 
ing to the wind and nuclear burning (IKato & Hachisu l l 1994). 

The decay phase of novae can be well represented 
with a sequence of ste ady state solutions as described by 
IKato & Hachisu I (119941) . We have solved a set of equations, 
i.e., the equations of motion, mass continuity, radiative dif- 
fusion, and conservation of energy, from the bottom of the 
hydrogen-rich envelope through the photosphere. The winds 
are accelerated deep inside the photosphere, so they are called 
"optically thick winds." 

2.2. The Reduced Opacity 

We assume that the opacity is effectively reduced by a fac- 
tor s, 

K eS = K/s, (1) 
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in the super-Eddington p hase, where k is the OPAL opacity 
(Iglesias & Rogers 1996) and s is the opacity reduction factor 
that represents the reduced ratio of the effective opacity in a 
p orous envelope. 

iKato & Hachisu I {2005) assumed that s is a function of the 
temperature and time, i.e., s is unity in the outer region of 
the envelope (log T < 4.7), but takes a certain constant value 
So (> 1) at the inner region logT > 5.0, and it changes lin- 
early between these values. Here sq is a function of time 
that has the maximum value at the optical peak and then 
gradually decreases to unity. Choosing an appropriate s{T,t), 
Kato & Hachisu reproduced the light curve of VI 974 Cyg 
in the super-Eddington phase. Once the temperature depen- 
dence of s is given, so is uniquely determined by fitting with 
both the optical and UV light curves. In the present paper, 
we first adopt the same temperature dependence of s as in 
IKato & Hachisul (120051) . We call it Model 1 . 

We adopt another type of function for s. We call it Model 2, 
in which s changes in a more inner region of the envelope, i.e., 
s = 1 at \ogT < 5.25 but s = so at \ogT > 5.45 and changes 
linearly between these values. 

The function of s should be closely linked with radiation in- 
stabilities against a porous structure of the atmosphere. How- 
ever, we do not know how and when the porous structure de- 
velops in a nova envelope and how much the opacity is re- 
duced. Therefore, in the present paper, we assume the two 
functions of s, i.e., Model 1 and Model 2. Model 1 corre- 
sponds to the case that the porous structure develops from the 
bottom of the envelope to a lower temperature region beyond 
the peak of the OPAL opacity at log J ~ 5.2. Therefore, the 
peak value of the opacity is reduced by a factor of so- Model 
2 corresponds to the case that the porous structure does not 
extend to the opacity peak at log T ~ 5.2. 

2.3. Optical Light Curve 

In VI 500 Cyg, which is one of the brightest novae, the 
temporal evolution of the spectrum and the fluxes are well 
understood as blackbody emission d uring the first three days 
and as free-free emission after that dGallagher & Nev1 ll976t 
lEnnis et al. lll977HDuerbeck & Seitter 111979b . In the model- 
ing of the super-Eddingt on phase we divide the optical light 
curve into three phases dKato & Hachisu Il2005l) . The first is 
the super-Eddington phase, in which we simply assume that 
photons are emitted at the photosphere as a blackbody with a 
photospheric temperature of T p h- In the next phase, the optical 
flux is dominated by free-free emission of the optically thin 
ejecta outside the photosphere. The flux of free-free emission 
can be roughly estimated as 

/r°° M 2 M 2 
NMdVtx / -aM^roc— sf-, (2) 

where F\ is the flux at the wavelength A, N e and Ni are the 
number densities of electrons and ions, V is the volume of 
the ejecta, M w i„d is the wind mass-loss rate, and v p h is the 
velocity at the photosphere. Here, we use the relation of 
Pwind =Mwind/47rr 2 v wind , and p wind and v wind are the density 
and velocity of the wind, respectively. We substitute M w ind, 
R p h and v p h from our best fit model. We cannot uniquely spec- 
ify the proportional constant in equation <(2J because radiative 
transfer is not calculated outside the photo sphere. Instead, 
we choose the constant to fit the light curve (lHachisu & Kato I 
l2(Ml2006h . 

When the nova enters a nebular phase, strong emission lines 
such as [O III] dominantly contribute to the visual light curve. 




JD 2,444,000+ (day) 

FIG. 1 . — UV light curve-fitting for V693 CrA. Calculated 1455 A con- 
tinuum fluxes are plotted against time for the WD models of Mwd =11 
Mq (dotted line), 1.2 M Q (dashed line), 1.25 M (thin solid fine), 1.3 M Q 
(thick solid line), and 1.35 Mq (dash-dotted line). The chemical composition 
of the envelope is assumed to be X — 0.35, Y — 0.33, Xcno = 0-2, Xn c =0.1, 
and Z = 0.02. The theoretical flux F1455 is calculated for an arbitrarily as- 
sumed distance of 1.0 kpc and no absorption (scale in the right-hand-side). 
Open circles denote the 1455 A continuum flux taken from Cassate lTa et al. I 
<200l . 

Then the visual light curve gradually deviates from our free- 
free light curve of equation (0. This is the third phase. 

2.4. UV1455A Light Curve 

After the optical maximum, the photospheric radius of the 
envelope gradually decreases while the photospheric temper- 
ature (Tph) increases with time. As the temperature increases, 
the main emitting wavelength of radiation shifts from op- 
tical to UV. The UV 1455 A flux reaches a maximum at 
log r p h ~ 4.4. After the UV flux decays, the supersoft X-ray 
flux finally increases. 

Figure Q] shows the theoretical 1455 A light-curves of var- 
ious WD masses with an envelope chemical composition of 
X = 0.35, Y = 0.33, Xqno = 0.2, X Ne = 0.1, and Z = 0.02. 
The evolutional timescale depends strongly on the WD mass. 
More massive WDs evolve faster than less massive WDs. 
The evolution speed also depends on the chemical composi- 
tion because enrichment of heavy elements drives more mas- 
sive winds through the opacity enhancement, which acceler- 
ate nova evolutions. Therefore, the duration of 1455 A burst 
depends on the WD mass and chemical composition of the 
envelope. 

Figure [2] depicts the dependence of such duration of the 
1455 A outburst for various WD masses and chemical com- 
positions. Here, the duration is defined by the full width at the 
half m aximum of the 1455 A light curve (lHachisu & Kato I 
2006). Once the chemical composition is determined, we can 
e stimate the WD mass fro m this figure. 

lHachisu & Kato I (|2006) fitted theoretical light curves with 
the 1455 A observation and determined the WD masses for 
VI 668 Cyg and VI 974 Cyg. They find that their WD masses 
show good agreement in the light-curve fittings of optical, in- 
frared, and X-ray for the entire period of the outburst. In this 
sense, the 1455 A light curve is a good indicator of the WD 
mass. In the present paper, we estimate the WD mass from 
the 1455 A light-curve fitting and use it in modeling of the 
super-Eddington phase. 

3. V693 CRA (NOVA CORONAE AUSTRALIS 1981) 

Nova V693 CrA wa s discovered by Honda (Kozai & Kosai I 
I198U ICaldwell|[T982h on 1981 April 2 near maximum at an 
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TABLE 1 
Model Parameters 





object 


VoVj CrA 




\M Q*7/1 /~\.<r 

v iy /4 t-yg 




V looo L.yg 


V J j 1 r up 


And 


outburst year 


1981 




1992 




1978 


1991 


1986 


opacity Model (1 or 2) 


1 


2 


1 


2 


2 


2 


2 


M W D (M Q ) 


1.3 


<_ 


1.05 


<_ 


0.95 


1.0 


1.0 


X 


0.35 


<- 


0.46 


<— 


0.45 


0.35 


0.45 


Y 


0.33 


<— 


0.32 


<— 


0.18 


0.23 


0.18 


XcNO 


0.2 


<— 


0.15 


<— 


0.35 


0.3 


0.35 




0.1 


<— 


0.05 


<— 





0.1 





Z 


0.02 




0.02 




0.02 


0.02 


0.02 


E(B-Vf 


0.2 




0.32 


«- 


0.40 


0.72 


0.25 


so at peak 


2.2 


2.7 


5.0 


7.2 


9.0 


5.5 


6.0 


Distance (kpc) 


4.4 


4.4 


1.8 


1.8 


4.1 


2.1 


4.3 


UV FWHM C (days) 


10 


10 


46 


46 


53 


54 


47 


Lmax (10 38 ergs" 1 ) 


2.68 


2.65 


4.4 


4.2 


4.0 


3.3 


3.2 


Mv.mtx (mag) 


-7.30 


-7.30 


-7.75 


-7.73 


-7.73 


-7.43 


-7.56 


mi/,max (mag) d 


6.5 


6.5 


4.6 


4.6 


6.6 


6.4 


6.4 


excess of super-Edd (mag) 


0.84 


0.84 


1.7 


1.7 


1.9 


1.5 


1.6 


duration of super-Edd (days) 


6 


6 


18 e 


I6 C ' 


16 


9 


12 


?n time (days) 


13 


13 


34 e 


32 e 


25 


23 


25 


AMqec (10- 5 M Q ) 


2.0 


2.7 


3.9 


4.8 


5.8 


3.1 


3.9 



a IVaru andingham et al. ( 1997) for V693CrA; Chochol et al.] fl997l) fo r V1974 Cyg; IHachisu & Katol <2006l) for 
V1668 Cyg; Saizar et al. 11996) for V351 Pup; Schwarz et al. 1 1997) for OS And b opacity reduction factor in 
equation {TJ at the optical maximum 1 -' calculated from our theoretical UV 1455 A light curve d apparent magnitude 
corresponds to M{J eak in the light curve fitting 11 estimated from the part of free-free light curve 



0.33, and X Ne = 0.20 dWilliams et al. II1985L taken from Table 
6 in Valandinghamet al. 1997), orX = 0.1 6, Y = 0.18, X CN o = 
0.36, and X Ne = 0.27 (I Andrea et al. 111994 . Considering these 
scattered values, we adopt X = 0.35, Y = 0.23, Xcno = 0.22, 
X^t e = 0.10 and Z = 0.02 in our model calculation. The theoret- 
ical light curve hardly changes if we increase the neon abun- 
dance from 0. 1 to 0.2 and decrease the helium from 0.23 to 
0.13. This is because the exchange of neon with helium does 
not affect either the hydrogen burning rate or the opacity. 

3.1. UV Light Curve and Distance 

FigureQ]depicts the IUE continuum UV fluxes obtained by 
ICassatella et al. I (120021) for the 1455 A band with a AA = 
20 A width (centered on A = 1455 A). The corresponding 
theoretical light curves are also plotted for five WD masses 
of 1.1, 1.2, 1.25, 1.3, and 1.35 M Q . Here, we assume the 
OPAL opacity (i.e., 5=1 throughout the envelope). For more 
massive WDs the evolution is faster and the UV flux decays 
more quickly. Both the 1.25 and 1.3 Mq WDs are consistent 
with the observation. We may exclude Mwd < 1-2M Q and 
Mwd > 1.35M . In the present paper,we adopt the 1.3 Mq 
m odel for later calculatio n. 

IVanla ndingha met al. I (Q997) obtained the reddening of 
V693 CrA to be E(B-V) = 0.2 ± 0. 1 mainly from the compar- 
ison with nova LMC 1990 No. 1 and also the comparison with 
the reddening of globular clusters within 10 degrees from the 
nova. They also suggested that the reddening is small because 
no interstellar absorption feature was seen in their spectrum. 
In the present paper, we assume E(B-V) = 0.2. 

Using this reddening, we estimate the distance to the star. 
The absorption at A = 1 455 A is ca lculated to be A\ = 
8.3 E(B-V) = 1.66 (e.g.. lSeatonlll979h . The observed peak 
flux is 7.64 x 10~ 12 erg cm" 2 s _1 A _1 whereas the theoretical 
peak value is 7.02 x 10" 10 erg cm" 2 s _1 A _1 for the 1.3 Mq star 
for a distance of 1 kpc. Therefore, the distance is calculated 
to be D = y/7.02 x 10" 10 /7.64 x lO" 12 /^ 1 - 66 / 2 - 5 ) = 4.4 kpc. 

In previous works, the distance to V693 CrA has been es- 
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FIG. 2.— Full width at the half maximum (FWHM) for the 1455 A light 
curve vs. the WD mass. Dash-dotted line: neon novae (X = 0.65, Y = 
0.27, Z C no = 0.03, X Ne = 0.03). Solid line: neon novae (X = 0.55, Y = 
0.30, X CN0 = 0.10, X Ne = 0.03). Dashed line: CO novae (X = 0.55, Y = 
0.23, Z C NO = 0.20). Solid line: CO novae (X = 0.45, 7 = 0.18, X CN o=0.35). 
Dotted line: CO novae (X = 0.35, Y = 0.33, X C no = 0.30), which is almost 
overlapped with the lower solid line. We assume Z = 0.02 for all the models. 
Arrows indicate the WD mass for each object. 

apparent magnitude of 6.5. The discovery magnitude was 
once reported to be 7.0 but wa s later correcte d to be 6.5 
(Kozai, private communication in Caldwell 1982). ISion et ail 
( 1986) suggested a white dwarf as massive as those in O-Ne- 
Mg novae from the broad emission line widths, high ejection 
velocities, large mass-loss rates, and th e presence of strong 
neon lines ( Vanlandingha met al. 1119971) . 

The abundance of V693 CrA was estimated from the IUE 
spectra to be X = 0.38, Y = 20, X cm = 0.15, and X Ne = 0.26 
dVanlandingham et al. I \l99% . X = 0.26, Y = 0.30, X CN0 = 
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FIG. 3.— Light-curve fitting for V693 CrA 1981. Thick solid lines: cal- 
culated V light curves for the blackbody assumption. Thin solid lines: calcu- 
lated 1455 A light curves. The WD mass is assumed to be 1.3 Mq with the 
envelope chemical composition of X = 0.35, Y = 0.33, Xcno = 0-2, Xj^ c = 0.1 
and Z = 0.02. Optical data are taken from AAV SO (crossesl. lCaldweH 1198 11 
squares), WE VFES tCassatella etaTll2004al. asterisks), and IAUC 3591, 
3590, 3594, and 3604 (small open circles and arrows for upper limit obser- 
vation before the outburst). The maximum magnitude on JD 2,444,697 was 
corrected from 7.0 to 6.5 I Caldwell 1982). The 1455 A data are the same as 
in Fig. [T] The opacity reduction factor .so is plotted by a thin solid curve in 
a linear scale between 1.0 (bottom) and 2.2 (at the peak) in the upper panel 
(Model 1), and 1.0 (at the bottom) and 2.7 (at the peak) in the lower panel 
(Model 2). The distance of 4.4 kpc is assumed, which is estimated from the 
1455 A light-curve fitting. 

timated to be as la rge as 8-12 kpc (ICaldwelll fT98ll 119821: 
lBroschlll981l, [l982). mainly because they assumed large ab- 
solute magnitudes of My = -8.75 to -10 from the absolute 
magnitude-rate of decline relations. This relation, however, 
is not very accurate for a single nova (Vanlandingha nTet al. I 
1 1997b and may overestimate the absolute magnitude. We will 
see that the peak magnitude is as faint as My = -7.30 in both 
Models 1 and 2 for the obtained distance of 4.4 kpc. We sum- 
marize our fitting results in Table Q] 

3.2. Optical Light Curve 

Figure[3]shows the theoretical light curves as well as the ob- 
servational data in optical. Because of few data points around 
the maximum and of scattered data after that, it is difficult to 
identify a shape of the optical light curve. Our best fit model 
for each opacity reduction factor shows good agreement with 
both the UV and optical light curves. These two models have 
similar properties as shown in Table 1 . 

In Figure [3] the thick solid curve denotes the visual light 
curve for the blackbody photosphere. The super-Eddington 
phase lasts 6 days in both models. Free-free emission gradu- 
ally becomes dominant as the photospheric temperature rises 
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FIG. 4.— Light-curve fitting for V1974 Cyg 1992 for Model 2. (a) Up- 
per panel. Thick solid line: V-magnitude from the blackbody photosphere. 
Dashed line: V-magnitude from the free-free emission calculated from equa- 
tion {2). Thin solid line: the opacity reduction factor sq in the linear scale 
between 1.0 (at the bottom) and 7.2 (at the peak). Optical data are taken 
from AAVSO (dots), ( b) Lowe r panel. Open circles: The 1455 A data taken 
from Cassatella et al. 1 2004b). In the upper panel the distance of 1.8 kpc is 
assumed, which is obtained from fitting in the lower panel. 

to logT > 4.0 and our theoretical curve for blackbody emis- 
sion deviates from the observation. Here, we do not plot the 
light curve of free-free emission phase, because we cannot 
exactly determine the constant in equation (f2]i for such large 
scattering data. 

In the case of V693 CrA the super-Eddington phase almost 
ends before the 1455 A flux rises. Therefore, our estimates 
of the WD mass and distance, which are determined from the 
UV flux fitting, are probably independent of our modeling of 
the super-Eddington phase. 

4. V1974 CYG (NOVA CYGNI 1992) 

iKato & Ha chisu (2005) presented a light curve model of 
the super-Eddington phase for a 1.05 Mq WD with a chem- 
ical composition of X = 0.46, Xqno = 0.15, X^ e = 0.05 and 
Z = 0.02. This is the prototype of Model 1. Here, we have 
calculated light curves for Model 2 with the same other pa- 
rameters as in Kato & Hachisu. Our best fit model in Figure|4] 
is very similar to that for Model 1 . This means that we cannot 
determine which opacity reduction factor s is preferable for 
VI 974 Cyg from the light-curve fitting. 

The distance is calculated from the UV light-curve fitting. 
We obtained 1.8 kpc with an extinction of E(B-V) = 0.32 
(IChochol et al.lfl997h . the same value as in Model 1. Our 
distance is sligh t ly larg er than the estimate of 1.7 kpc by 
lHachisu & Kato I (120051) obtained from the UV light-curve fit- 
ting using the normal opacity (s = 1) models, because the lu- 
minosity at the UV peak is still super-Eddington (sq = 1.15 in 
Model 1 and s Q = 1 .2 in Model 2 at the UV peak), and the 
UV flux is 15-20% larger than that in the normal opacity 
models. Our distance is consistent with those discussed in 
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FIG. 5.— Same as FigH but for V1668 Cyg 1978. (a) Upper panel. 
Thick solid line: V-magnitude from the blackbody photosphere for Model 
2. Thin solid line: the opacity reduction factor sq in the linear scale be- 
tween 1.0 (at the bottom) and 9.0 (at the peak). The WD mass is assumed 
to be 0.95 Mq with X = 0.45, X C NO = 0.35, and Z = 0.02. Op tical data 
are taken from AAVSO (dots) and Mallama & Skillman] 119791 crosses). 
Data of the Stri i mgren y band magnitude (crosses with a circle) are from 
Gallagh er et al. 1 098 01. (b) Lower panel. UV data (open circles) are taken 
fromlCassatella et alTl 420021) . The distance of 4.1 kpc is assumed in the up- 
per panel, which is obtained from the 1455 A light-curve fitting in the lower 
panel. 

IChochol et al. I (119971) . a most probable value of 1.8 kpc. 

5. V1668 CYG (NOVA CYGNI 1978) 

VI 668 Cyg wa s discovered on 1978 September 10.24 UT 
dMorrisonlll978l) . two days before its optical maximum of 
mv,max = 6 .04. This object was als o well observed with IUE 
satellite. lHachisu & Kato I (120061) presented a light curve 
model from shortly after the optical peak until the end of the 
outburst. Their best fit model is a WD mass of 0.95 M Q with a 
chemical composition of X = 0.45, Xcno = 0.35, and Z = 0.02. 
In the present paper, we adopt their parameters and reproduce 
the super- Eddington phase. 

VI 668 Cyg shows a similar but slightly steeper light curve 
around the peak compared with V1974 Cyg. As shown in 
Figures [4] and [5] there is a remarkable difference in the 1455 
A fluxes. In V1668 Cyg the 1455 A flux remains low in 
several days before rising up, whereas in VI 974 Cyg it once 
decreases before the optical maximum and rises again. This 
means that the photospheric temperature in VI 668 Cyg re- 
mains low in a relatively long time around the optical peak. 
This difference led us to find no best fit solutions for Model 1 . 
Figure|5]shows our best fit light curve for Model 2, where the 
WD mass a nd the chemical compos ition is assumed to be the 
same as in Hachisu & Kato I d2006l) . This model represents 
well both the UV and y magnitude light curves. 

5.1. Distance from UV Light-Curve Fitting 

Figure [5] shows that the 1455 A flux decays quickly around 
JD 2,443,820, which can be attributed to the formation of an 
optically thin dust shell. iGehrz et al. I (H980) reported an ex- 
cess of the infrared flux that peaks at JD 2,443,815 whereas 
no significant drop in the visual magnitude is observed. Even 



if the dust shell is optically thin, the grain condensation leads 
to a redistribution of UV flux into infrared. Considering this 
effect, which is not included in our theory, our 0.95 M Q WD 
model shows good agreement with the 1455 A observation. 

It is difficult to estimate the interstellar absorption in the 
direction of V1668 Cyg becau se there are only a few star 
well photometrically observed. Hachi su & Kato I (120061) re- 
examined the distance-reddening law in the direction of 
VI 668 Cyg and estimated the redden ing to be E(B-V) = 0.4 
and obtained the distance of 3.6 kpc. Stickl and et al. I (119811) 
estimated E(B-V ) = 0.4 ± 0. 1 from 2200 A feature. We adopt 
E(B-V) = 0.4. 

The distance is derived to be 4.1 kpc with E(B-V) = 0.4. 
This value is somewha t larger than 3.6 kpc estimated by 
lHachisu & Kato I ([2006) with the normal opacity. In our 
model, the opacity reduction factor is still larger than unity 
(so = 1.45 at the 1455 A peak) as shown in the upper panel 
of Figure|5] Both the bolometric an d UV fluxes are enhance d 
compared with those in the model o flHachisu & Katol (120061) . 
T herefore, we obtain a la rger distance. 

ISlovak & Vogt I (fl9 79) estimated the distance to be 3.3 kpc. 
iDuerbeck et al. I (119801) obtained a much smaller distance of 
d = 2.3 kpc based on the same stars. It is very difficult to es- 
timate the distance because of the small number of stars and 
the very patchy Ay (r) re l ation in the direct i on of t he nova (see 
lHachisu & Kato1l2006l) . iGallagher et all (119801) derived the 
distance t o be 4.4 kpc adopting M y = -7.8 for E(B-V) = 0.3. 

119801) preferred a smaller distance 



However, Galla gher et al. 
of d ~ 2 kpc, which was derived assuming that the peak lu- 
minosity is equal to the Eddington luminosity because a star 
with the super-Eddington luminosity may be dynamically un- 
stable. This argument cannot be applied to our model, how- 
ever, because the envelope is settled down into a steady-state 
even in the super-Eddington phase. Stickl arid et al. I d 198 lb 
suggested the distance to be 2.2 kpc by equating the max- 
imum luminosity and the Eddington luminosity. They also 
obtained a distance of 3.6 kpc from the relation between the 
maximum magnitude and the rate of decline. They did not 
take this larger distance because the large acceleration of mat- 
ter cannot be expected in the optically thin region. However, 
this is not the case of our optically thick wind. Considering 
these arguments, we conclude that our distance of 4. 1 kpc is 
reasonable. 

Using the distance of 4.1 kpc obtained from the 1455 
A light-curve fitting, we derive the peak magnitude of 
My.max = —7.73, i.e., super-Eddington by 1.9 mag. 

5.2. Optical Light Curve 

Our theoretical V-magnitudes are shown in the upper panel 
of Figure|5] This curve is placed so as to satisfy the distance 
modulus of (m-M) v = 1.24 + 51og(4.1kpc/10pc)= 14.30. 

There are rich observational data of visual magnitude. 
How ever, we focus on the Stromgren y band light curve (taken 
from lGallagher et al. II 19801 denoted by the crosses with a cir- 
cle), which lies along the bottom edge of the visual data, be- 
cause the y filter is designed to avoid strong emission lines 
in the nebular phase and reasonably represents the continuum 
fluxes of novae. We regard that our light curve follows the y 
magnitude because our model ignores such line contributions. 

The nov a enters a coronal phase 53 days after the optical 
maximum dKlare et aill 19801) and th e spectrum shows strong 
nebular emission lines dKaler|[l980l) . After that, the differ- 
ence between visual and y magnitudes becomes significant as 
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shown in Figure [5] 

The spectral development of VI 668 Cyg can be understood 
with our model as follows. The spectrum near the optical 
maximum was reported to be similar to that of an interme- 
diate F star (lOrto lani et al.l l 19781) . and also be consistent with 
the principal spectrum which is characterized by weak hydro- 
gen emission li nes and absorption l ines of neutral and singly 
ionized metals (Kla re & Wolf1ll978l) . In our model, the pho- 
tospheric temperature is as low as log T p ^ = 3.87, most of hy- 
drogen is recombined in the region around the photosphere. 
Therefore, these properties are consistent with the observed 
features. When the magnitude declined by abo ut 0.8 mag 
(three days after the maximum), iKlare et aT] (1 19801) reported 
that the nova shows the diffuse enhanced spectrum in which 
strong emission features dominate. In our model the photo- 
spheric temperature rises to log T p ^ = 3.92 at this time and the 
ionization degree of hydrogen is quickly increasing with time. 
This is consistent with the appearance of strong H/3 emission. 
The no va enters a corona l phase 53 days after the optical max- 
imum (IKlare et ai.lll980h . In our model the temperature rises 
to log r p h = 4.57 at this time and the photon flux is dominated 
by UV. 

The characteristic properties of the light curve are summa- 
rizes in Table 1. The mass of the envel ope expelled in the 
outburst is estimate d to be 2 x 10~ 5 M^ jGehrz et al. I [l 9801) 
and 5.5 x 10" 5 M Q (IStickland et al. 11198 lh . Our model gives 
an ejecta mass of 5.8 x 10~ 5 M Q , which is roughly consistent 
with the observational estimates. 

6. V351 PUP (NOVA PUPPIS 1991) 

V351 Pup was discovered on 1991 December 27 by 
Camilliiri] (fl992h near maximum. The light curve shown 
in Figure [6] resembles that of VI 668 Cyg in Figure [5] In 
our modeling we adopt X = 0.35, Y = 0.23, Xqno = 0.30, 
X Ne = 0.10 and Z = 0.02 after Saizar et al.'s (1996) estimate 
from the WE spectra, i.e., X = 0.36, Y = 0.24, X CN o = 0.26, 
and X^ e = 0.12. The reddening is obta ined to be E(B- V) = 
0.79-0.92 from emission-line ratios (Wil liams! 1 19941) and 
E(B - V) = 0.72 ± 0.1 from ratios of recombination lines 
(ISaizar et al. 1119961) . Here we take E(B-V) = 0.72. 

As optical data are poor around the maximum and no y 
magnitude data were reported, we assume Model 2 as in 
V1668 Cyg and a slightly larger WD mass of 1.0M Q . The 
resultant light curve is shown in Figure [6] which shows good 
agreement with both the 1455 A and optical light curves. The 
distance is estimated to be 2.1 kpc from the UV light-curve 
fitting. 

7. OS AND (NOVA ANDROMEDAE 1986) 

OS And w as discovered by Suzuki on 1986 December 5 
dKosai 111986b . The optical light curve in Figure|7]shows a 1.5 
mag dip that lasts about 30 days owing to dust formation. Cor- 
responding to this dip, th e 1455 A light curve sh ows a quick 
decrease at JD 2,446,800 (ICassatella et al. 120021) . Apart from 
the deep dip, the optical and 1455 A light curves resemble to 
those of V1668 Cyg and V351 Pup in the first 20 days. More- 
over, the UV spectr um is very similar to that of V351 Pup in 
the first two weeks (Sonnebor n et al.l["l992l) . 

With such a strongly absorbed UV light curve, it is difficult 
to search for a best fit model. Therefore, we assume the same 
WD mass and the same s function as those in V351 Pup. The 
chemical composition is assumed to be the same as in VI 668 
Cyg, regarding OS And as a CO nova. Note that the difference 
in the composition does not make large difference as shown 
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FIG. 6.— Same as Figg] but for V351 Pup 1991. (a) Upper panel. Thick 
solid line: V-magnitude from the blackbody photosphere. Thin solid line: 
the opacity reduction factor so in the lineal' scale between 1.0 (at the bottom) 
and 6.0 (at the peak). Op tical data are taken from AAVSO (dots), IUE VFES 
Cassate llTetaU J2004a. asterisks), and IAUC 5422, 5423, 5427, 5430, 5437, 
5447, 5455, 5493, 5503, and 5527 (squares and arrows), (b) Lower panel. 
Data of UV 1455 A continuum are taken from lCassatella et al. I <2002l) . The 
distance of 2. 1 kpc is assumed in the upper panel, which is obtained from the 
1455 A light-curve fitting in the lower panel. 
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FIG. 7.— Same as Figg] but for OS And 1986. (a) Upper panel. Thick 
solid line: V-magnitude from the blackbody photosphere. Thin solid line: 
the opacity reduction factor .so in the linear scale between 1.0 (at the bot- 
tom) and 6.0 (at the peak). Optical data are taken from AAVSO (dots), 
Kikuch ilFaTl (1981 squares), and IAUC 4281, 4282, 4286, 4293, 4298, 
4306, 4342 and 4360 (asterisks), (b) Lower panel. The 1455 A data are 
taken from Cassatella et alT] 420021) . The distance of 4.3 kpc is assumed in 
the upper panel, which is obtained from the fitting in the lower panel. 
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FIG. 8. — The wind mass-loss rate CM w j nt i), envelope mass (AM), and the 
ratio of the gas pressure to the total pressure at the critical point of steady-state 
wind (/3cr) are plotted against the opacity reduction factor, so, in the super 
Eddington phase. The right-end point of each curve corresponds to the optical 
peak while the left-end point corresponds to sq = 1, i.e., the epoch when the 
super Eddington phase ends. Time goes on from right to left. Dashed/Solid 
line with squares: Model 1/Model 2 of V693 CrA, respectively. Both of them 
are almost overlapped in all three cases of M w i n d, AM, and f) a , Dashed/Thin 
solid line: Model 1/Model 2 of V1974 Cyg. Thick solid line: V1668 Cyg. 
Dash-dotted line: V351 Pup. Dotted line: OS And. 

in the two light curves of VI 668 Cyg and V351 Pup. The 
blackbody light curve in Figure|7]shows good agreement with 
both the visual and UV data in the first 20 days. 

The distance is estimated to be 4.3 kpc from the 1455 
A light curve fi tting in the lower pa nel for a reddening of 
E(B-V) = 0.25 (ISchwarz et al. Ill997h . Our value of 4.3 kpc 
is roughly consistent with 5.1 ± 1.5 kpc obtained from a com- 
parison of UV fl uxes between OS And and Nova LMC 1992 
(ISchwarz et aT1l997h . 

8. DISCUSSION AND SUMMARY 

Table 1 summarizes the model parameters and our main 
results, i.e., from top to bottom row, (1) object name, (2) 
outburst year, (3) adopted function of the opacity reduction 
factor, (4) WD mass, (5)-(9) adopted chemical composition, 
(10) reddening, (11) opacity reduction factor so at the opti- 
cal maximum, (12) distance estimated from the peak of the 
1455 A light-curve, (13) duration of the 1455 A outburst de- 
fined by the full width at the half maximum (FWHM), (14) 
peak bolometric luminosity, (15) absolute and (16) apparent 
V-magnitudes corresponding to the peak luminosity, (17) ex- 
cess of the super-Eddington in V-magnitude, i.e., the differ- 
ence between the peak magnitude of our model (My max ) and 
the peak magnitude of a light curve model with the normal 
opacity (s = 1) for the same WD mass and the same envelope 
composition, (18) duration of the super-Eddington phase, (19) 
time in which My drops by 3 magnitude from the peak in our 
theoretical model, (20) mass ejected during the period from 
the first point of each light curve until the wind stops. 

8.1. Opacity Reduction Factor 



There are no time-dependent calculations of nova outbursts 
in which porous instabilities widely develop. So, it is difficult 
for us to estimate how much the opacity is reduced in the nova 
envelopes. Here, we have simply assumed s to be a function 
of the temperature and time. Then, we have determined so so 
as to reproduce the observed light curve of each object. 

The referee kindly pointed out that the opacity reduction 
factor would be expressed as a function of the current state of 
the envelope which is independent of time. After a porous in- 
stability widely develops in the envelope, it may settle into an 
equilibrium state. If our opacity reduction factor represents 
such an envelope state, it may be a function of a small num- 
ber of physical parameters/variables that represent the current 
state of the envelope. 

In order to search for such parameters/variables, we plot the 
wind mass-loss rate (M w indX envelope mass (AM), and ratio 
of the gas pressure to the total pressure at the critical point of 
each steady wind solution (/?„) against So in Figure[8] We plot 
seven light-curve models for five objects in Table 1. These 
seven models are different from each other in their character- 
istic properties such as the WD mass, chemical composition, 
and opacity model. 

First of all, we can easily see that the difference between 
Model 1 and Model 2 hardly makes large difference on these 
three variables; the curves for Model 1 and Model 2 are very 
close to each other in V693 CrA and VI 974 Cyg. 

Next, we see that the two models of V693 CrA are separated 
from the others in the wind mass-loss rate. This is because 
the WD of V693 CrA is as massive as 1.3 M Q whereas the 
other WDs are 0.95- 1.05 Mq. The wind mass-loss rate is 
much larger on a massive WD than on a less massive one if 
we compare them a t the same envelope m ass, as shown in 
our previous work dKato & Hachisulll994h . Therefore, the 
difference in the wind mass-loss rates of V693 CrA in Figure 
[8]can be attributed to the difference in the WD masses. So we 
may conclude that the mass-loss rate itself is not a main factor 
that determines the opacity reduction factor so. 

In the envelope mass, all the curves are almost similar to 
each other. This suggests that s is closely related to the en- 
velope mass. In other words, the porous instability involves 
at least a wide area of the envelope, and so decreases as the 
envelope mass decreases with time due to wind mass-loss. 

In the ratio of the gas pressure to the total pressure at the 
critical poin t of steady-state wind so lutions, where the wind is 
accelerated (Kato & Hachisu 1994), we see that all the curves 
are also very similar to each other. The ratio [3 is almost con- 
stant in a region below the critical point. Therefore, (3 a - rep- 
resents a mean value of /? in the region where the opacity is 
reduced. 

We have checked other physical variables such as the lu- 
minosity, photospheric radius, temperature, and radius at the 
critical point, but we found that all these variables are largely 
scattered from each other, i.e., not bunched like in AM and 
f3 a - of Figure [8] 

We may conclude that the envelope mass and (3 cr are the 
key parameters that represent the current state of envelopes 
with a porous structure. Probably, /?„ should be closely linked 
with the property that the porous instability is a kind of radi- 
ation instability and the envelope mass may be related with 
the property that the porous instability involves a large part 
of the envelope, although we do still not know the accurate 
condition for the porous instability. 
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8.2. White Dwarf Mass and Chemical Composition 

Both the WD mass and chemical composition are fixed in 
the present work. If we choose a different set of these param- 
eters, we have different time scales of theoretical light curves 
not only in the very early phases but also in the later phases 
of the outbursts including the epochs when the wind mass- 
loss stops and hydrogen burning ends. For V1974 Cyg and 
VI 668 Cyg, multiwavelength observations are available un- 
til the very late phase of the outburst and their WD masses 
are determined consistently wi th these observations (see, e.g., 
lHachisu & Katol l200l 12006. for VI 974 Cyg and VI 668 
Cyg). Therefore, we adopt their estimates. 

The decline rates of nova light-curve depend strongly on the 
WD mass and weakly on X, very weakly on Xcno but hardly 
on X^ e . Therefore, the largest ambiguity in the estimation of 
WD masses comes from the accuracy of hydrogen content X. 
The dependency of the WD mass on X is roughly estimated 
as 

Mwd<7) «M W d(0.55) + 0.5(X-0.55), (3) 

when 0.35 < X < 0.65 and 0.03 < X CN o < 0-35, using the 
"univ ersal decline law" of nova light curves (lHachisu & Kato I 
2006) that nova light curves are almost homologous except 
for the very early phase (i.e., the super-Eddington phase dis- 
cussed here). Here, Mwd(0.55) means the WD mass esti- 
mated for X = 0.55 and the original value of X = 0.55 can 
be replaced with any other value, for example, X = 0.35. 

For the other three objects, i.e., V693 CrA, V351 Pup, 
and OS And, there is no observational data in the late phase, 
i.e., when the wind stopped and when the hydrogen burning 
ended. Then we determined the WD mass only from the UV 
light curve fitting with a fixed chemical composition as shown 
in Figure [T] If we adopt a different set of the chemical com- 
position, we have a slightly different WD mass as can be esti- 
mated from equation (0 or Figure [2] 

8.3. Distance 

Nova distances are determined from the comparison be- 
tween observed UV fluxes and calculated fluxes at the UV 
peak. In case of V693 CrA and V351 Pup, the super- 
Eddington phase ended, i.e., sq = 1, at the UV peak as shown 



in Figures [3] and [6] So, the distance can be determined in- 
dependently of the reduced factor of s. In case of V1974 
Cyg, the super-Eddington phase still continues at the UV 
peak. The distances derived for Model 1 and Model 2 are 
the same within the accuracy of two digits as shown in Table 
1, although these values themselves are somewhat larger than 
the values derived without the super-Eddington phase (s = 1) 
as already mentioned in §5.1. From these results, it may be 
concluded that the opacity reduction factor does not affect so 
much the distance estimate. This is because, in all the five 
objects, their super-Eddington phases had already or almost 
ended (so — 1 0) at the UV peak. 

8.4. Summary 

Our main results are summarized as follows; 

1. We present light-curve models of the super-Eddington 
phases for five IUE classical novae based on the optically 
thick wind theory, with an a ssumption tha t the opacity is re- 
duced in a porous envelope (Shaviv 2002). Our models rea- 
sonably reproduce the optical and 1455 A light curves. 

2. The duration of the 1455 A light curve is a useful indi- 
cator of the WD mass, especially when the chemical compo- 
sition is known. 

3. The distance is derived from the comparison between the 
observed peak value of the 1455 A flux and the corresponding 
calculated value unless a dust shell absorbs the UV flux at its 
maximum. The derived distances of the five IUE novae are 
consistent with the previous estimates. 
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